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Abstract 

Starting from the presence of non-ideal plasma effects due to strongly 
coupled plasma in the early stage of relativistic heavy-ion collisions, 
we investigate, from a phenomenological point of view, the relevance 
of non-conventional statistical mechanics effects on the rapidity spec- 
tra of net proton yield at AGS, SPS and RHIC. We show that the 
broad rapidity shape measured at RHIC can be very well reproduced 
in the framework of a non-linear relativistic Fokker-Planck equation 
which incorporates non-extensive statistics and anomalous diffusion. 



1 Introduction 

It is expected that hadrons dissociate into a plasma of their elementary con- 
stituents, quarks and gluons (QGP), at a density several times the nuclear 
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matter density and/or at temperatures of the order of T c = 170 MeV, which 
is the critical temperature of the transition from the hadronic gas phase 
to the QGP phase, as expected from lattice QCD calculations. In central 
Au+Au collisions at RHIC energy densities are reached that are far above 
the critical energy density (of the order of 1 GeV/fm 3 ) for a transition to the 
QGP phase. 

Since interactions among quarks and gluons become weak at small dis- 
tance or high energy, one usually expects that QGP is a weakly interacting, 
ideal plasma, which can be described by perturbative QCD. However, this is 
rigorously true only at very high temperature, while non-perturbative phe- 
nomena prevail up to temperatures of several times T c PQ. For example, as 
it is known from several approaches, partons in a hot plasma acquire a tem- 
perature dependent Debye mass, which deeply modifies the QCD dynamics, 
yet inducing collective effects. Hence the final states in an ultrarelativistic 
nucleus-nucleus collision can be very sensibly affected by non-ideal plasma 
effects, including long-range interactions and memory effects. Many proper- 
ties concerning the formation of the expected new phase and the consequent 
hadronization process are still under debate (21 El EJ E] • Pre-thermalization or 
metastability conditions can be realized and the standard equilibrium statis- 
tical mechanics assumptions can not be taken for granted in the description 
of the system toward equilibrium jB] . 

In this paper, we want to investigate the relevance of such non-ideal and 
non-equilibrium plasma effects within a phenomenological study of the net- 
baryon rapidity distribution in various relativistic heavy ion collisions, up to 
the central Au+Au collisions at the highest RHIC energy of ^5^=200 GeV. 
The interest of this observable lies in the fact that in these processes nuclear 
matter reaches high energy densities and nuclei undergoing central collisions 
strongly reduce their original longitudinal momentum. This loss of rapidity, 
usually referred to as baryon stopping, is an important characteristic to un- 
derstand the reaction mechanism at high energy density and the net rapidity 
distribution is very sensitive to the dynamical and statistical properties of 
nucleus-nucleus at high energy. 

We study the evolution of the rapidity distribution from a macroscopic 
point of view by using a non-linear relativistic Fokker-Planck equation and 
we will show that the observed broad rapidity shape could be a signal of non- 
equilibrium properties of the system in which non-extensive Tsallis statistical 
mechanics |H] emerges in a natural way. Similar approaches have been 
considered in the past to analyze transverse momentum distributions [2*ll§l llU[ 
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El El and power law spectra at large p± in terms of various non conventional 
extensions of the Boltzmann-Gibbs thermostatistics, including the Tsallis 
one. 

The net-proton rapidity distribution appears to be a more elusive observ- 
able with respect to equilibrium considerations, since the main component is 
expected to be related to the initial beam protons; however the results ob- 
tained at RHIC do not appear easily reproduced within standard approaches. 
Encouraging results have been obtained by Wolschin |T3J El [E] within a 
three-components Relativistic Diffusion Model (RDM). Here again a possi- 
ble link with Tsallis statistics is suggested [T4"| . 

The paper is structured as follows. In Sec. 2, we qualitatively review 
some basic features of a non-ideal QGP like the one hypothetically gener- 
ated in the relativistic heavy-ion collisions at AGS, SPS and RHIC energies. 
We also show that the main assumptions contained in the derivation of the 
standard dynamical kinetic equations, describing the evolution of the system 
toward the equilibrium, are no longer valid. In Sec. 3, we introduce a non- 
linear relativistic kinetic equation containing anomalous diffusion effects in 
the framework of the non-extensive Tsallis thermostatistics. In Sec. 4, we 
study the net proton rapidity spectra comparing the AGS, SPS and RHIC 
data. Conclusions are reported in Sec. 5. 



2 Non-ideal QGP plasmas and kinetic assump- 
tions towards the equilibrium 

In the literature, an ordinary plasma is usually characterized by the value of 
the plasma parameter T ^B] 

r = M (i) 

defined as the ratio between potential energy (U) versus kinetic energy (T). 
When r C 1, one has a dilute weakly interacting gas; the Debye screening 
length Ad is much greater than the average interparticle distance r and a 
large number of particles is contained in the Debye sphere. Binary colli- 
sions induced by screened forces produce, in the classical case, the standard 
Maxwell-Boltzmann velocity distribution. If T « 0.1 4- 1, then Ad ~ tq, 
and it is not possible to clearly separate individual and collective degrees 
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of freedom: this situation refers to a weakly interacting, non-ideal plasma. 
Finally, if T > 1, the plasma is strongly interacting, Coulomb interaction and 
quantum effects dominate and determine the structure of the system. 

The quark-gluon plasma close to the critical temperature is a strongly in- 
teracting system. In fact, following Ref . [3 El H] , the color-Coulomb coupling 
parameter of the QGP is defined, in analogy with the one of the classical 
plasma, as 

r « c-^- , (2) 

47rr T ' V ; 

where C = 4/3 or 3 is the Casimir invariant for the quarks or gluons, respec- 
tively; for typical temperatures attained in relativistic heavy ion collisions, 
T ~ 200 MeV, a s = # 2 /(4vr) = 0.2 0.5, and r ~ n' 1 ^ ~ 0.5 fm (n being 
the particle density for an ideal gas of 2 quark flavors in QGP). Consequently, 
one obtains T ~ 1.5 — 5 and the plasma can be considered to be in a non-ideal 
liquid phase El • 

In these conditions, the generated QGP does not satisfy anymore the 
basic assumptions (BBGKY hierarchy) of a kinetic equation (Boltzmann or 
Fokker-Planck equation) which describes a system toward the equilibrium. 
In fact, near the phase transition the interaction range is much larger than 
the Debye screening length and a small number of partons is contained in the 
Debye sphere E] • Therefore, the collision time is not much smaller than 
the mean time between collisions and the interaction is not local. The binary 
collisions approximation is not satisfied, memory effects and long-range color 
interactions give rise to the presence of non-Markovian processes in the ki- 
netic equation, thus affecting the thermalization process toward equilibrium 
as well as the standard equilibrium distribution. In the next section we will 
see how such effects can be taken into account, from a macroscopic point of 
view, by means of an appropriate generalization of the standard statistical 
mechanics. 



3 Relativistic non-extensive kinetic equations 

In many-body long-range-interacting systems, it has been recently observed 
the emergence of long-living quasi stationary (metastable) states character- 
ized by non-Gaussian power law velocity distributions, before the Boltzmann- 
Gibbs equilibrium is attained. In this case, the standard statistical mechanics 
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is no longer appropriate to describe such a behavior [T7| [TBI HH] . Indeed, the 
basic assumption of the Boltzmann-Gibbs statistical mechanics is that the 
system can be subdivided into a set of non-overlapping subsystems, the total 
entropy of which is the sum of the entropies of the independent subsystems. 
In the presence of memory effects and long range forces the entropy, which 
is a measure of the information about the particle distribution in the states 
available to the system, does not satisfy the extensivity property. 

Recently, there is an increasing evidence that the generalized non-extensive 
statistical mechanics, proposed by Tsallis [HE], can be considered as an ap- 
propriate basis for a theoretical framework to describe physical phenomena 
where long-range interactions, long-range microscopic memories and/or frac- 
tal space-time constraints are present. A considerable variety of physical 
issues show a quantitative agreement between experimental data and the- 
oretical models based on Tsallis' thermostatistics. In particular, there is a 
growing interest to high energy physics applications of non-extensive statis- 
tics. Several authors outline the possibility that experimental observations 
in relativistic heavy-ion collisions can reflect non-extensive statistical me- 
chanics effects during the early stage of the collisions and the thermalization 
evolution of the system GJJ Q2 ES3 I2H 122 I2H| • 

In order to study from a phenomenological point of view experimental 
observables in relativistic heavy-ion collisions, we can introduce the basic 
macroscopic variables in the language of relativistic kinetic theory following 
the Tsallis' prescriptions for the non-extensive statistical mechanics. In this 
framework the particle (four-vector) flow can be generalized as [21] 



where we have set % — c — 1, x = x^ — (t, x), p = p^ = (p°,p) and 
p° = a/p 2 + m 2 is the relativistic energy. In the above Z q = J dfl [f(x,p)] q is 
the non-extensive partition function, <El stands for the corresponding phase 
space volume element and q is the deformation parameter. The limit q — > 1 
corresponds to ordinary statistics. The four- vector iV M = (n, j) contains the 
probability density n = n(x) (which is normalized to unity) and the proba- 
bility flow j = j(x). The energy- momentum tensor contains the normalized 





and the energy-momentum flow as 
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g-mean expectation value 1 of the energy density, as well as the energy flow, 
the momentum and the momentum flow per particle. 

On the basis of the above definitions, one can show that it is possible to 
obtain a generalized non-linear relativistic Boltzmann equation [21] 

p^[f(x,p)] g = C q (x,p), (6) 

where the function C q (x,p) implicitly defines a generalized non-extensive 
collision term 

C 9 (x,p) = l0^ d ^{h q [f\f\]W(p',p\\p, Pl ) 
zj p l p p l 

-h q [f,fi}W(p,p 1 \p',p' 1 )} . (7) 

Here W(p,pi\p' ,p\) is the transition rate between a two-particle state with 
initial four-momenta p and p\ and a final state with four-momenta p' and 
p\; h q [f, fi] is the g-correlation function relative to two particles in the same 
space-time position but with different four-momenta p and pi, respectively. 
Such a transport equation conserves the probability normalization (number 
of particles) and is consistent with the energy-momentum conservation laws. 
The collision term contains a generalized expression of the molecular chaos 
and for q > implies the validity of a generalized if-theorem, if the following, 
non-extensive, local four-density entropy is assumed 

d 3 p 

where we have used the definition ln g x = (x 1-9 — l)/(l — q). At equilibrium, 
the solution of the above Boltzmann equation is a relativistic Tsallis-like 
(power law) distribution and can be written as 

1/(1-9) 



W = ~K [ ^-P»fi(x,p)n\n q f(x,p) - 1] , (8) 



f €q (p) = TT 
Z 9 



1 - (1 - q) 



pHJ^ 
k B T 



(9) 



In the limit q — > 1, Eq.Q reduces to the standard relativistic equilibrium 
Juttner distribution. 

The above relations represent the basic framework in which, in the next 
section, will be studied the net-baryon rapidity distribution near equilibrium. 

1 The g-mean expectation value is denned as [HI I24| : 

(0(x)>, = / d 3 pO(x,p)[f(x,p)}«/ [ dn[f(x, P )}" . (5) 
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4 Net-baryon rapidity distribution 



The energy loss of colliding nuclei is a fundamental quantity in order to 
determine the energy available for particle production in heavy-ion collisions. 
Since the baryon number is conserved and rapidity distributions are only 
affected by rescattering in the late stages of the collision, the measured net- 
baryon (B — B) distribution retains information about the energy loss and 
allows one to obtain the degree of nuclear stopping. 

Recent results for net-proton rapidity spectra in central Au+Au collisions 
at RHIC [22] show an unexpectedly large rapidity density at midrapidity in 
comparison with analogous spectra at lower energy at SPS j2H| and AGS |27j . 
As outlined from different authors, such spectra can reflect non-equilibrium 
effects even if the energy dependence of the rapidity spectra is not very well 
understood [23I2EI- 

We are going to study the evolution of the rapidity distribution from 
a macroscopic point of view by means of a non-linear relativistic Fokker- 
Planck equation which can be viewed as the near-equilibrium approximation 
of the generalized non-linear Boltzmann equation (JSj). Therefore, the use 
of the Fokker-Planck equation appears to be correct only not too far from 
equilibrium (13 HH 113 1221 EOl Ell E21 E3 EH 

In order to study the rapidity spectra, it is convenient to separate the ki- 
netic variables into their transverse and longitudinal components, the latter 
being related to the rapidity y. If we assume that the particle distribution 



function f(y } m_|_, t), at fixed transverse mass m± = ym 2 + p]_, is not appre- 
ciably influenced by the transverse dynamics (which is considered in thermal 
equilibrium), the non-linear Fokker-Planck equation in the rapidity space y 
can be written as 



d r „, , n d 

di [f{y ^ t)] = dy- 



J{y, m ± ) [f(y, m ± , t)} + D-^-[f(y, m ± , *)]" 



(10) 



where D and J are the diffusion and drift coefficients, respectively. 

Tsallis and Bukman [201 have shown that, for linear drift, the time depen- 
dent solution of the above equation is a Tsallis (non-relativistic) distribution 
with /i = 2 — q and that a value of q ^ 1 implies anomalous diffusion, i.e., 
\y{t) — yhiif)} 2 scales like t a , with a = 2/(3 — q). For q < 1, the above equa- 
tion implies anomalous subdiffusion, while for q > 1, we have a superdiffusion 
process in the rapidity space. 
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Let us note that a similar approach, within a linear and non-linear Fokker- 
Planck equation, has been previously studied in Ref. [T3J CH HH] involving 
directly the time evolution of the rapidity distribution instead of the parti- 
cle distribution function f(y,m±,t), as in the Eq. ffTUj) . The two approaches 
are equivalent only if the particle distribution is completely decoupled in 
the transverse and in the longitudinal coordinates and the drift coefficient 
J does not depend on transverse momentum. This is the case, for example, 
when the drift coefficient is assumed to be linear in the rapidity space. In 
Ref. fHH ITU I15j. the author first study a linear Fokker-Planck equation with 
a linear drift coefficient and a free parameter diffusion coefficient (connected 
to a collective longitudinal expansion) which implies a strong violation of the 
fluctuation-dissipation theorem. This approach fails to reproduce the rapid- 
ity spectra at RHIC energy unless it is assumed that a fraction of net protons, 
near midrapidity, undergoes a fast transition to local thermal equilibrium. 
The non-linear Fokker-Planck equation for the rapidity distribution func- 
tion is also considered in Ref. [T51 ITi} ITS] , even if only approximate solutions 
are considered by means of a linear superposition of power-law distribution 
functions. 

We claim that the choice of the diffusion and the drift coefficients plays 
a crucial role in the solution of the above non-linear Fokker-Planck equa- 
tion (fTU|) . Such a choice influences the time evolution of the system and its 
equilibrium distribution. This important point has been extensively studied 
in Ref.jS] where the authors obtained a generalized fluctuation-dissipation 
theorem, which implies an implicit relation between the drift and diffusion 
coefficient in the framework of a linear non-relativistic Fokker-Planck equa- 
tion. 

We are going to show that by generalizing the Brownian motion in a 
relativistic framework, the standard Einstein relation is satisfied and Tsal- 
lis non-extensive statistics emerges in a natural way from the non-linearity 
of the Fokker-Planck equation. In fact, by imposing the validity of the Ein- 
stein relation for Brownian particles and setting henceforward the Boltzmann 
constant k B to unity, we can generalize to the relativistic case the standard 
expressions of diffusion and drift coefficients as follows 

J D = 7 T, J(y,m ± ) = 7 m ± smh(y) = 7 p\\ , (11) 

where pii is the longitudinal momentum, T is the temperature and 7 is a 
common constant. Let us remark that the above definition of the diffusion 
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and drift coefficients, previously introduced by us in Ref.[20j, appears as the 
natural generalization to the relativistic Brownian case in the rapidity space. 
The drift coefficient which is, as usual, linear in the longitudinal momentum 
p\\ becomes non-linear in the rapidity coordinate. 

It is easy to see that the above coefficients give us the Boltzmann sta- 
tionary distribution in the linear case (g = fi — 1), while the equilibrium 
solution f eg (y,m±) of Eq. lfTUJl . with \i = 2 — q, is a Tsallis-like distribution 
(introduced in Eq.(|9*|l) with the relativistic energy E = m± cosh(y) 



Out of equilibrium the rapidity distribution at fixed time can be obtained by 
means of numerical integration of Eq.(|l(J|) with delta function initial condi- 
tions depending upon the value of the experimental projectile rapidities. The 
rapidity distribution at fixed time is then obtained by numerical integration 
over the transverse mass m± as follows 



where m is the mass of the considered particles and c is the normalization 
constant, fixed by the experimental data. The rapidity spectra calculated 
from Eq. (J13|) will ultimately depend on two parameters: the "interaction" 
time r int = yt and the non-extensive parameter q. 

Let us observe that in the numerical solution of Eq.(13) we have not 
explicitly assumed the presence of longitudinal flow. This fact does not 
exclude the dynamical effects of a collective (mainly due to rescattering) 
flow but rather incorporates a description of it in the adopted non-extensive 
statistical mechanics. In fact, as explicitly shown in Ref. [2] by studying the 
transverse mass spectrum, dynamical collective interactions are intrinsically 
involved in the generalized statistical mechanics and, in a purely thermal 
source, a generalized g-blue shift factor (depending on q and on m±) appears. 
We expect a similar behavior also in the longitudinal degrees of freedom of 
the system. Let us also notice that collective transverse flow effects in the 
framework of a non-extensive statistical mechanics have been investigated in 
Ref. [Hj as well. 

In Figs. 1 and 2 we plot the numerical solution of Eq.s (fTU|) and (JT3*j) at 
different interaction times Tj n4 in the linear case (q = 1) and in the non-linear 
case, with q = 1.25, keeping the same initial conditions. In comparing Fig. 1 




1/(1-9) 



(12) 




(13) 
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with Fig. 2, we can observe that the rapidity spectra appear to be broader in 
the non-linear case. This is a consequence of the anomalous (super) diffusion 
process implied in the time evolution of Eq. (|10j) and in the nature of the 
power law distribution function which, for q > 1, enhances the probability 
to find particles at high rapidity values. 
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Figure 1: Numerical solution of Eq.s ljlO)) and (fTHj) at different interaction 
times Ti n t in the linear case (q — 1). 



In Fig. 3, we report the obtained rapidity distribution (full line) for the 
net proton production (p — p) compared with the experimental data of RHIC 
(Au+Au at = 200 GeV, j2SJ), SPS (Pb+Pb at ^/JJ^ = 17.3 GeV, 

|26j) and AGS (Au+Au at ^snn = 5 GeV, [21] )• The parameters employed 
for the three curves are: q = 1.485 with Ti nt = 0.47 for RHIC, q = 1.235 
with Ti n t = 0.84 for SPS and q = 1.09 with T int = 0.95 for AGS, respectively. 
We notice that, although q and T int appear, in principle, as independent 
parameters, in fitting the data they are not. Indeed, we can see that only in 
the non-linear case (q ^ 1) there exists one and only one (finite) time T int for 
which the obtained rapidity spectrum well reproduces the broad experimental 
shape. On the contrary, for q = 1, no value of T int can be found, which allows 
to reproduce the data. 

We obtain a remarkable agreement with the experimental data by increas- 
ing the value of the non-linear deformation parameter q as the beam energy 
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Figure 2: The same of Fig.l but in the non-linear case with q = 1.25. 



increases. At AGS energy, the non-extensive statistical effects are negligible 
and the spectrum is well reproduced within the standard quasi-equilibrium 
linear approach. At SPS energy, non-equilibrium effects and non-linear evo- 
lution become remarkable (q = 1.235) and such effects are even more evident 
for the very broad RHIC spectra (q = 1.485). Let us observe that such an 
excellent agreement with the RHIC experimental data has not been reached 
in the similar (but different, as pointed out in the introduction) approach of 
Ref. [T3*l CHI especially in reproducing the experimental points far from 
the midrapidity region. 

From a phenomenological point of view, we can read the larger value of 
the parameter q for the RHIC data, corresponding to non-linear anomalous 
(super) diffusion, as a signal of the non-ideal nature of the plasma formed at a 
temperature larger than the critical one. As confirmed by recent microscopic 
calculations [3J E] , strongly coupled non-ideal plasma is generated at energy 
densities corresponding to the order of the critical phase transition temper- 
ature and in such a regime we find, in our macroscopic approach, strong 
deviations from the standard thermostatistics. At much higher energy, such 
as LHC, we can expect a minor relevance of such non-ideal effects since the 
considerable energy density reached is far above the critical one. However, 
an anomalous diffusion behavior should still affect the time evolution process 
toward the equilibrium due to long range color magnetic forces which remain 



11 



unscreened (in leading order) at all temperatures. 
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Figure 3: Rapidity spectra for net proton production (p—p) at RHIC (Au+Au 
at = 200 GeV, BRAHMS data), SPS (Pb+Pb at ^IJ^ = 17.3 GeV, 

NA49 data) and AGS (Au+Au at ^/IJ^ = 5 GeV, E802, E877, E917). 



5 Conclusions 

We have studied the rapidity spectra from a macroscopic point of view 
by means of a non-linear kinetic equation which preserves the fluctuation- 
dissipation theorem by introducing an appropriate relativistic generalization 
of the drift and diffusion coefficients. Such a generalized evolution equation 
lies inside the framework of Tsallis' non-extensive thermostatistics and con- 
tains anomalous diffusion effects which are strongly related to the presence 
of non-Markovian memory interactions and long-range color forces. Multi- 
particle rescattering, very large at SPS and at RHIC, could be a signal of 
this effects. We recall that non-extensive features in relativistic collisions are 
suggested in several works, even if a microscopic justification of these effects 
is still lacking and lies out of the scope of this paper. The relevant point 
which was found in the present work is that the larger the beam energy of 
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the heavy ions, the greater are the values of the deformation parameter q, 
which are required to reproduce the experimental rapidity spectra. This be- 
havior of the g-parameter can be viewed as a phenomenological indication of 
the onset of a strongly coupled QGP plasma phase in the early stages of the 
collision. 
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